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Abstract 

A new phenomenological technique for using constant amplitude loading data to predict fatigue life from a 
variable amplitude strain history is presented. A critical feature of this reversal-by-reversal model is that the 
damage accumulation is inherently non-linear. The damage for a reversal in the variable amplitude loading 
history is predicted by approximating that the accumulated damage comes from a constant amplitude loading 
that has the strain range of the particular variable amplitude reversal. A key feature of this approach is that 
overloads at the beginning of the strain history have a more substantial impact on the total lifetime than 
overloads applied toward the end of the cycle life. This technique effectively incorporates the strain history in 
the damage prediction and has the advantage over other methods in that there are no fitting parameters that 
require substantial experimental data. The model presented here is validated using experimental variable 
amplitude fatigue data for three different metals. 
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1. Introduction 

One of the most prevalent methods for testing 
the fatigue properties of a material is to construct 
a constant amplitude strain-life curve. In this 
method, a sample is strained cyclically between two 
strain levels until failure. This test is performed on 
a number of identical samples at different strain 
magnitudes, and the results are plotted as a strain- 
life curve. This - 27V/" curve can be fitted to 
the Basquin-Manson-Coffin, Q, i, Q (BMC) equa- 
tion, 
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where ^ is half of the strain range, Ae e is the 
clastic strain, Ae p is the plastic strain, a / is the 
fatigue strength coefficient, E is the elastic mod- 
ulus, b is the fatigue strength exponent, e / is the 
fatigue ductility coefficient, c is the fatigue ductil- 
ity exponent, and 27V/ is the number of reversals 
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until failure. The low cycle fatigue regime is con- 
sidered a product of plastic deformation, and the 
high cycle portion is related primarily to clastic de- 
formation. This approach gives a reasonably accu- 
rate prediction of sample life at constant amplitude 
cyclic strains. 

However, parts in service are rarely subjected 
to idealized constant amplitude cyclic strains and 
instead undergo variable amplitude loading. The 
Palmgren-Miner (PM) rule is a widely used ap- 
proach for predicting part lifetime under variable 
amplitude loading. This method hypothesizes that 
the damage caused by each stress state in a variable 
amplitude load history is a function of the number 
of times that the particular cyclic stress state oc- 
curs and the number of cycles it would take for the 
sample to fail from a constant amplitude history at 
that stress state. Q Mathematically the Palmgren- 
Miner rule is written, 
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where the sum is taken over all stress states with 
i denoting a particular state, rii is the number of 
cycles at the ith stress state, Nfi is the number of 
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cycles to failure if the sample is cycled under con- 
stant amplitude loading at the ith stress state, and 
C is a constant. Based on the assumptions made 
by the model, C should be 1; experimentally it is 
found to range between 0.7 and 2.2. This variability 
is evidence of the failure of the Palmgren-Miner rule 
to accurately predict fatigue lifetimes. Although it 
is well known that the Palmgren-Miner rule is in- 
accurate, its conceptual simplicity and the minimal 
amount of data necessary for implementation makes 
it a popular method for estimating fatigue life. 

Improving upon the Palmgren-Miner rule has 
been a major focus of researchers studying variable 
amplitude fatigue. A good review of cumulative 
damage and life prediction theories through the end 
of the last century is presented by Fatemi and Yang. 

In recent years models have been developed 
based on a variety of techniques that couple the- 
ories of fracture mechanics and empirical observa- 
tions. 0, 0, 0, M, [H H El, E3, El Other meth- 
ods take into account the residual stresses caused 
by the plasticity of the material at the crack tip, 
and crack tip closing phenomena, 17, 18 1 Al- 



to i = 2N T , 
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though these modeling methods are more accurate 
than the Palmgren-Miner rule, they require sub- 
stantially more experimental data to fit the neces- 
sary parameters. 

In this paper a new method, free of fitting pa- 
rameters, is demonstrated for estimating strain life 
under variable amplitude loading. It is unique from 
other models in that the only data used for input 
is the constant amplitude strain-life curve and the 
cyclic stress-strain response. Using this model it 
is possible to accurately predict the variable am- 
plitude strain life of specimens using a relatively 
small amount of experimental data that can eas- 
ily be generated. In the section following this in- 
troduction the analytical details of the model are 
presented. Next the experimental procedure and 
model implementation are explained. In the results 
and discussion section the measured and predicted 
strain lifetimes are presented and compared. The 
theoretical lifetimes predicted from this model are 
compared to lifetime predictions presented in the 
literature. In the final section the paper is suc- 
cinctly summarized. 



2. Model Details 

In this model the total damage, Dt, is the sum of 
the damage of all reversals, Z),, ranging from i = 1 
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In this definition Di is the normalized damage 
caused by the ith reversal and failure occurs when 
Dt = 1. The damage caused by each reversal is 
determined using the well-known constant ampli- 
tude strain-life relation. Using this approach, the 
strain history is incorporated as damage accumu- 
lation. The damage accrued in each step is calcu- 
lated using a relatively simple algorithm and con- 
stant amplitude strain-life data. 

Following examples from Ref. the damage 
during fatigue is assumed to be due to a single crit- 
ical crack propagating across the width of the spec- 
imen. A good description of the crack growth per 
reversal, for a constant amplitude strain, is given 
by the hyperbolic sine function. For a crack size, a, 
the rate of crack tip advance, after 2N k reversals, 
can be written 



da . (2N k 
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where 2Nf is the total number of reversals to fail- 
ure for the given strain amplitude, and p is a scal- 
ing factor, which will be discussed in more detail 
later. Eq. (QJ is a natural expression of damage per 
reversal, and has the same functional form as has 
been used to describe changes in crack growth rates 
previously. It accurately reflects the phenom- 

ena in that during the initial stages of damage the 
rate of crack propagation is low and as damage is 
accrued the rate of crack growth increases. Other 
expressions for the rate of crack tip advance, which 
arc phcnomcnologically similar, are presented else- 
where and will not be discussed here. [ij], [H| 

Using Eq. (j4]) for damage, the normalized damage 
due to the ith reversal is expressed, 



Di 



E^ 1 +1 sinh(jfp)-E^sinh(if 

Ej2i sinh 
JX T+I sinh(f^) dN k 



f" f sinh (ffip) dNj 

(5) 

where Nt is the number of reversals required to 
achieve the accumulated damage, Dt, were the 
damage due to a constant amplitude strain range, 
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Ap- • The number of reversals to failure at this con- 
stant strain range is 2Nf. The denominator nor- 
malizes the damage per reversal such that at failure 
the total damage is 1. 

The current state of damage, Dt, is known and 
is expressed as the sum of incremental damage, Di, 
from Eq. ©, 
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where Nt is the number of damaging reversals to 
cause the damage Dt at the constant amplitude 
strain 4/. The summations in Eqs. and ^ can 
be approximated as integrals, which allows Eq. ([5]) 
to be solved to find 



2N T = I —f- ) x 
cosh -1 
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The general approach for calculating a specimen 
life under variable amplitude loading goes as fol- 
lows. Begin by determining the initial total dam- 
age. This is typically near zero, if the part begins 
pristine and the strain-life curve is well defined at 
all strains of interest. A case for non-zero starting 
damage will be discussed later in section |U The 
strain range for the first tensile reversal is 4p , and 
D\ is calculated using Eq. ([5]) and the BMC re- 
lation, Eq. ([T]), which is fit to constant amplitude 
strain life data. The value of D\ is added to Dt- 
The next tensile strain is The lifetime at this 
strain range along with Dt, are used to calculate 
Nt from Eq. ([7]), which is used in Eq. ([5]) to cal- 
culate Z?2, which is again added to Dt- This pro- 
cess continues until Dt = 1, at which point failure 
is predicted. An algorithm for implementing this 
model is demonstrated in section |U 



3. Experimental Data 

Published experimental data are used to validate 



steel and compared their experimental results to an 
effective strain damage model based on the work of 
DuQuesnay. [2(| Colin and Fatemi published ex- 
perimental data for 304L stainless steel and 7075 
T6 aluminum. [2lj 

The load histories used in the calculations were 
recreated from the descriptions detailed by Pereira 
ct al., and Colin and Fatemi. [HHH The P355NL1 
specimens were subjected to a variety of load 
blocks, including a high to low scheme, a low to high 
scheme, a low to high to low scheme, and random 
loading, examples of which are shown in Figs. [TJ [3J 
[3J and|U The 304L stainless steel fatigue-life results 
were prepared for periodic, fully-reversed overloads, 
shown in Fig. [5j and random loading. The 7075 T6 
aluminum samples were subjected to random load- 
ing. Unlike the shaped loading blocks, the random 
load history was not reproduced from the litera- 
ture, but instead a strain history file was created 
using a random number generator, filtering the ran- 
dom number stream to ensure that each iteration 
reversed the strain. A representative sample of the 
random reversal data is shown in Fig. 01 

H-L load block (100 cycles) 




the model. [20|, |2l| Pereira et al. tested P355NL1 



100 

N, reversals 

Figure 1: An example of a high to low strain loading 
block after the work of Pereira et al. 1201 



4. Model Implementation 

The model described in section [2~] above was im- 
plemented using the algorithm detailed in this sec- 
tion and applied to the strain history data discussed 
in section[31 To begin, the first two strains from the 
strain history were converted to stresses using the 
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L-H load block (100 cycles) 



Example random load block 




100 150 
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Figure 2: An example of a low to high strain loading 
block after the work of Pereira et al. |2(| 



L-H-L load block (100 cycles) 




Figure 3: An example of a low to high to low strain 
loading block after the work of Pereira et al. i2Qj 



Ramberg-Osgood stress-strain relationship, 
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where K and n' were fit from the cyclic stress- 
strain behavior. The Morrow mean stress correc- 
tion Q, 



(9) 



where a m = a <-+ 1 + CJ > i s the mean stress, was used to 
calculate Nf. The value of Nf was used with Dt 
to calculate Nt from Eq. (0. Finally, the damage 
from this reversal, Di, was determined from Eq. (J5j) . 



Figure 4: An example of random loading. Note 
that the random number stream was filtered to en- 
sure that each iteration reversed the direction of the 
applied strain. 

Periodic overload (fully reversed) 




20 30 

N, reversals 

Figure 5: An example of a fully reversed periodic 
overload loading block after the work of Colin and 
Fatemi. EOl 



The damage was added to the total damage, Dt- 
If the total damage was greater than 1, then the 
specimen was deemed to have failed due to this re- 
versal, otherwise the procedure was continued using 
the new value of Dt and the next strain taken from 
the strain history. This was repeated until failure 
occurs, when Dt > 1. 

The constant p, in Eq. ((4]), appropriately scales 
the incremental damage. For this implementation 

it was selected to be p = log (ZN]^ , where 2N f was 
from the BMC equation and 7 was taken as — c/2ej. 
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In this way, p scales with the applied strains and 
the incremental damage, da/dN, has the correct 
functional relationship to the strain amplitude. 

The integrated da/dN curve, which expresses to- 
tal damage as a function of the number of cycles, 
has a general shape that is well known, [jj @, 0] 
Careful inspection of experimental data of dam- 
age as a function of number of cycles reveals that 
a smooth well fitting curve does not always inter- 
cept the damage axis at 0. In many experiments, 
damage has been observed to accumulate rapidly to 
around 5 to 10% early in the specimen's life, and 
then slow to the da/dN crack propagation model 
that is well known. [lji @, @] To account for the 
rapid damage accumulation that occurs during the 
initial cycling, the starting damage was assigned to 
be Dt = 0.05 for all of the variable amplitude life- 
time prediction data presented here. 

5. Results and Discussion 

A comparison of the published experimental data 
and the predictions of this work can be seen in 
Figs. ISITlTfl In Figs. [5] and [7] the experimental data 
for P355NL1 steel, from Pereira et al., is shown for 
a maximum strain of 1.05% and 2.10% for the load- 
ing blocks discussed in section|3]above. [2(| The ex- 
perimental results are compared to those predicted 
from the present model in addition to the model 
of Pereira and DuQuesnay. The current model is 
in good agreement with both the experimental and 
theoretical results from Peirera et al. 

The experimental results from Colin and Fatcmi 
for 304L stainless steel are shown in Figs. [5] and 
[5] for the loading blocks discussed above at var- 
ious strain ranges. [2lj The results for 7075 T6 
aluminum are shown in Fig. 1101 The accuracy of 
the present model is compared to the predictions 
of the Palmgren-Miner rule both with and with- 
out the Smith- Watson- Topper correction. The cur- 
rent model is again in good agreement with the 
experimental results. It is worth noting that the 
model presented here agrees with the experimen- 
tal data both when the Palmgren-Miner results are 
non-conservative by an order of 2 or 3 and when 
they agree well with experiment. This is evidence 
that the results form the current work are more 
profound than a simple lifetime reduction from the 
Palmgren-Miner rule. 

The strength of the present model is due to the 
natural inclusion of the strain history when deter- 
mining the inflicted damage caused by a strain re- 
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Figure 6: A comparison of this work and the work of 
Pereira et al. for P355NL1 steel specimens tested by 
Pereira et al. and presented in Ref . [20| . A variety 
of load schemes are used at a maximum strain of 
1.05%. 
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Figure 7: A comparison of this work and the work of 
Pereira et al. for P355NL1 steel specimens tested by 
Pereira et al. and presented in Ref. [20| ■ A variety 
of load schemes are used at a maximum strain of 
2.1%. 



vcrsal. Both the effect of the immediately preceding 
strain and the effect of the relative age of a spec- 
imen are included. To calculate As a strain and 
the immediately preceding strain must be known 
to determine the applied strain range and the mean 
stress correction. The damage inflicted by a partic- 
ular strain reversal depends not only on the am- 
plitude, but also on the total state of damage at 
the instant of the reversal. Take the example of a 
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Figure 8: A comparison of the predictions of this 
work to the Palmgren-Miner rule (PM) with and 
without a Smith- Watson- Topper correction (SWT) 
prepared by Colin and Fatemi for 304L stainless 
steel specimens. 
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The data here represent spec- 
imens subjected to a constant amplitude fully re- 
versed strain history punctuated by periodic fully 
reversed overloads. 
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Figure 9: A comparison of the predictions of this 
work to the Palmgren-Miner rule (PM) with and 
without a Smith- Watson- Topper correction (SWT) 
prepared by Colin and Fatemi for 304L stainless 
steel specimens. [2l|The data here represent speci- 
mens subjected to random loading. 



periodic, fully-reversed, overload. If one considers 
the integrated da/dN curve for a constant ampli- 
tude cyclic strain, applying a fully reversed over- 
load would advance the position of the subsequent 
constant amplitude reversals on this curve substan- 
tially. From the shape of the curve, it is apparent 
that an overload early in the specimen's life that in- 
crements the damage will have a more substantial 
impact on the specimen's total life, compared to an 



Figure 10: A comparison of the predictions of this 
work to the Palmgren-Miner rule (PM) with and 
without a Smith- Watson- Topper correction (SWT) 
prepared by Colin and Fatemi for 7075-T6 alu- 
minum specimens. 2l| The data here represent 



specimens subjected to random loading. 



overload applied later. 

The more recent models, such as those of [3, OH 
[lH , include residual strains and crack-tip plasticity 
or crack tip closure phenomena. These models ac- 
count for the cumulative damage through localized 
plasticity near the crack tip, or a change in the ef- 
fective stresses due to crack tip closure. They are 
more accurate than the Palmgren-Miner model and 
other simple models because they account for the 
ordering of the applied strains and localized dam- 
age near the crack tip; however, they require sub- 
stantial experimental support. By comparison, the 
present model only includes materials information 
from the cyclic stress-strain curve and the constant 
amplitude strain life curve. It should be noted that 
for some analysis involving extreme load cases, a 
more sophisticated model than the one presented 
here may be necessary. In particular, cases where 
residual plasticity will play a dominant role in crack 
growth rates through strain hardening and residual 
stresses will likely require a model that explicitly 
deals with plasticity near the crack tip. 

6. Summary 

One of the greatest engineering challenges of the 
last 100 years is that of predicting the strain-life 
relationships of mechanical components undergo- 
ing variable amplitude loading. In spite of exten- 
sive studies, no conclusive model has been deter- 
mined. Although many useful models have been 
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developed, many require cumbersome amounts of 
experimental data. Here we report a new method, 
free of fitting parameters, for making an accurate 
variable amplitude strain-life prediction using ba- 
sic constant amplitude fatigue data. The method is 
validated using data from experiments performed 
on P355NL1 steel, 304L stainless steel, and 7075 
T6 aluminum. [2^, 21 The present model fits the 



experimental data well for a variety of load spec- 
tra and materials and the algorithm is simple to 
implement. 
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